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Abstract. We establish new interaction estimates for a system introduced 
by Baiti and Jenssen. These estimates are pivotal to the analysis of the wave 
front-tracking approximation. In a companion paper we use them to con¬ 
struct a counter-example which shows that Schaeffer’s Regularity Theorem for 
scalar conservation laws does not extend to systems. The counter-example we 
construct shows, furthermore, that a wave-pattern containing infinitely many 
shocks can be robust with respect to perturbations of the initial data. The 
proof of the interaction estimates is based on the explicit computation of the 
wave fan curves and on a perturbation argument. 


1. Introduction 

We deal with the system of conservation laws 

(1) dtU + d^[F^{U)\ = Q. 

The unknown U = U{t,x) attains values in 

U : [0,+oo[xIR —> 

/u 

{t,x) i-A U = \ V 

\w 

—> is defined as 

4:[{v - l)u - w] +ripi{U) \ 

\^4|u(u — 2)u — {v — l)w| + r]p 3 {U) j 

In the previous expression, the parameter r\ attains values in the interval [0, l/4[ 
and to simplify the exposition we fix the functions pi and pa by setting 


(3) 

Pi{U) = 2uw — 2u^{v 

(4) 

PsiU) = w‘^ — u^(v — 
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and the flux function E,, : 

/ 

(2) F,{U) := 


1 
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Note, however, that for some of the results discussed in the following the precise 
expression of the functions pi and pa is irrelevant. 

System (1),(2) was introduced by Baiti and Jenssen in [3, 19] and it was used 
to construct an example of a Cauchy problem where the initial data have finite, 
but large, total variation and the L°°-norm of the admissible solution blows up in 
finite time. More recently, the authors of the present paper used the Baiti-Jenssen 
system (1) to exhibit an explicit counter-example which shows that Schaeffer’s reg¬ 
ularity result for scalar conservation laws does not extend to systems, see [11]. The 
counter-example we construct shows, furthermore, that a wave-pattern containing 
infinitely many shocks can be robust with respect to perturbations of the initial 
data. We refer to § 2.1 in the present paper for a brief overview of these counter¬ 
examples. See also [10]. 

This note aims at establishing new quantitative interaction estimates for the 
Baiti-Jenssen systems (1),(2). The estimates we obtain are pivotal to the analysis 
of the so-called wave front-tracking approximation of the Cauchy problem obtained 
by coupling (1) with an initial datum 17(0,-) = Uq. We refer to [5, 14, 18] for 
an extended discussion on the wave front-tracking approximation. Here we only 
mention that the wave front-tracking algorithm is based on the construction of a 
piecewise constant approximation of the Cauchy problem. Under suitable conditions 
on the initial datum Uq and on the flux function F^, one can show that the wave 
front-tracking approximation converges to an admissible solution of the Cauchy 
problem, see in particular the analysis in [5]. In [11] we construct wave front¬ 
tracking approximations of the Cauchy problems obtained by coupling (1) with 
suitable initial data. We then rely on the wave front-tracking approximation to 
establish qualitative properties of the limit solutions. In the following we do not 
consider all the possible interactions one has to handle when constructing the wave 
front-tracking approximation. We only discuss those that we encounter in [11] and 
that cannot be handled by relying on straightforward considerations on the structure 
of the flux Fff. 

Before going into the technical details, we make some further remarks. First, in 
the present note we fix a very specific system in the wider class considered in [3]. 
The motivation for this choice is twofold: i) it simplifies the notation and ii) the 
analysis in the present note is sufficient for the applications in [11]. Note, moreover, 
that in the proof of Lemma 1.1 we use (although not in an essential way) the exact 
expression of the function evaluated at ij = 0. However, we are confident that 
our results can be extended to wider classes of systems of the type considered in [3] . 

Second, in this note the only occurrences where we explicitly use the precise 
expression of the functions pi and p^ is in the results discussed in § 2.3. More 
precisely, the proofs of Lemmas 1.1 and 1.2 both rely on a perturbation argument: 
we first show that our system verifies the statement of the lemmas in the case rj = 0 
and then we show that the same holds provided rj is sufficiently small. The proof of 
Lemma 1.2 is completely independent of the specihc expression of pi and p^. In the 
perturbation argument in the proof of Lemma 1.1 we use some results from § 2.3, 
but we never directly use the specific expression of pi and p 2 ■ 

Third, the Baiti-Jenssen (1) system in not physical, in the sense that it does not 
admit strictly convex entropies, see [3] for a proof. It is natural to wonder whether or 
not the results established in the present note can be extended to physical systems. 
Very loosely speaking, by combining Lemmas 1.2 and 1.1 below with the analysis 
in [11, §3.1-3.2] we get the following statement. Under suitable conditions, the only 
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waves generated at the interactions between two shocks are shock waves, or, more 
precisely, no rarefaction waves are generated at the interaction between two shocks. 
There are actually several physical systems that share this property: for instance, 
one can consider the 2x2 example discussed by DiPerna in [15, §5] and assume 
that the data have sufficiently small total variation. We refer to [ 6 , §4] for the 
analysis of shock interactions for this system. On the other hand, a much more 
challenging question is whether or not there is any physical system that exhibit 
the same behaviors as those discussed in [3, 11]. In other words, one can wonder 
whether or not a physical system can i) exhibit hnite time blow up or ii) violate 
the regularity prescribed, for scalar conservation laws, by Schaeffer’s Theorem. To 
the best of the authors’ knowledge, the answers to the above questions is presently 
open. 

We now give some technical details about the estimates we establish. First, we 
point out that the Baiti-Jenssen system (1) is strictly hyperbolic in the unit ball, 
which amounts to say that the Jacobian matrix DF^j admits three real and distinct 
eigenvalues 


(5) Ai(17) < A 2 (t/) < A3(17) 

for every U such that \U\ < 1. Also, if r? > 0 every characteristic field is genuinely 
nonlinear. In other words, let rq,..., ra denote the right smooth eigenvectors asso¬ 
ciated to the eigenvalues Ai, A 2 , A 3 . Then 

( 6 ) VA([/)-r*(t7)>c>0 


for some suitable constant c > 0 and for every i = 1,2,3 and \U\ < 1. In the 
following, we distinguish three families of shocks: we term a given shock 1 -, 2 - or 
3-shock depending on whether the speed of the shock is close to Ai, A 2 or A 3 . 

We also point out that establishing interaction estimates for system (1) boils 
down to the following. Consider the so-called Riemann problem, namely the Cauchy 
problem obtained by coupling ( 1 ) with an initial datum in the form 


(7) 


17(0, cc) 


Ui X < 0 
Ur X > 0, 


where Ue, Ur G are constant states. The above problem admits, in general, in¬ 
finitely many distributional solutions: we term admissible the solution constructed 
by Lax in the pioneering work [21], see § 2.2 for a brief overview. Establishing inter¬ 
action estimates for ( 1 ) amounts to establish estimates on the admissible solution of 
the Riemann problem (l)-(7) in the case when Ui and Ur satisfy suitable structural 
assumptions. 

The first case we consider is the case of the interaction of two 2-shocks, see 
Figure 1, left part. In other words, we assume that there is a state Um G such 
that 


• Ui and Um are the left and the right states of a Lax admissible 2-shock, 

• Um and Ur are the left and the right states of a Lax admissible 2-shock and 

• the shock between Ui and Um has higher speed than the shock between Um 
and Ur- 

We now give an heuristic formulation of our interaction estimate and we refer to § 3 
for the rigorous statement, which requires some technical notation. Here we only 
point out that the strength of a shock is a quantity defined in § 2.2 which is propor¬ 
tional to the modulus of the difference between the left and the right state of the 
shock. 
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Figure 1. Left: interaction between two 2-shocks. Right: inter¬ 
action between a 2-shock and a 1-shock. 




Lemma 1.1. Fix a constant a such that 0 < a < 1/2 and set := (a,0,—a). 
Consider the interaction between two 2-shocks and assume that the states Uf, and Ur 
are sufficiently close to C/jj. If the strengths of the interacting 2-shocks are sufficiently 
small, then the admissible solution of the Riemann problem (l)-(7) is obtained by 
patching together a 1-shock, a 2-shock and a 3-shock. 

We remark that the relevant point in the above result is that the solution of the 
Riemann problem that we consider in the statement contains no rarefaction wave. 

The second case we consider is the case of the interaction between a 1-shock and 
a 2-shock, see Figure 1, right part. In other words, we assume that there is a state 
Um G such that 

• Uf, and Um are the left and the right states of a Lax admissible 2-shock, 

• Um and Ur are the left and the right states of a Lax admissible 1-shock. 
The case of the interaction of a 3-shock with a 2-shock is analogous. We now give 
an heuristic formulation of our result and we refer to § 4 for the rigorous statement. 

Lemma 1.2. Consider the interaction between a 1-shock and a 2-shock and assume 
both shocks have suffieiently small strength. Then the admissible solution of the 
Riemann problem (l)-(7) is obtained by patching together a 1-shock, a 2-shock and 
a 3-shock. Also, we establish quantitative bounds from above and from below on the 
strength of the outgoing shocks, see formulas (35). 

Note that the fact that the three outgoing waves are shocks follows from the 
analysis in [3]. Also, the bound from above on the strength of the outgoing 3- 
shocks follows from by now classical interaction estimates, see [5, Page 133, (7.31)]: 
the main novelty in Lemma 1.2 is that we have a new bound from below on the 
strength of the outgoing 3-shock, see the left hand side of formula (35). This 
estimate is important for the analysis in [11]. 

This note is organized as follows. In § 2 we go over some previous results. 
In particular, in § 2.1 we provide some motivation for studying the Baiti-Jenssen 
system (1) by describing two counter-examples that use it. In § 2.2 we recall some 
results from [21] and in § 2.3 we apply these results to the Baiti-Jenssen system. 
In § 3 we discuss the interaction of two 2-shocks and in § 4 the interaction of a 
I-shock and a 2-shock. 
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2. Overview of previous results 

For the reader’s convenience, in this section we go over some previous results. 
More precisely: 

§ 2.1: we discuss two counter-examples based on the Baiti-Jenssen system (1): the 
original one in [3] and a more recent one devised in [11]. 

§ 2.2: we follow the famous work by Lax [21] and we outline the construction of 
the solution of the Riemann problem. 

§ 2.3: we apply Lax’s construction to the Baiti-Jenssen system. 

2.1. Counter-examples based on the Baiti-Jenssen system. This paragraph 
is organized as follows: 

§ 2.1.1: we discuss the counter-example in [3] 

§ 2 . 1 . 2 : we discuss the counter-example in [ 11 ]. 

Before dealing with the specific examples, we recall two main features of the Baiti- 
Jenssen system: first, it is strictly hyperbolic, namely (5) holds. Note that strict 
hyperbolicity is a standard hypothesis for results concerning systems of conservation 
laws, see [14]. Also, if 77 > 0 every characteristic field is genuinely nonlinear, which 
means that ( 6 ) is satisfied for every j = 1, 2, 3. This is a remarkable property because 
loosely speaking systems where all the characteristic field are genuinely nonlinear 
are usually better behaved than general systems. For instance, the celebrated decay 
estimate by Oleinik [23], which applies to scalar conservation laws with convex 
fluxes, has been extended to systems of conservation laws where all the characteristic 
field are genuinely nonlinear, see for instance the works by Glimm and Lax [17], by 
Liu [22] and, more recently, by Bressan and Colombo [7], Bressan and Goatin [ 8 ] 
and Bressan and Yang [9], while for balance laws we refer to Christoforou and 
Trivisa [12]. 

2.1.1. Finite time blow up of admissible solutions with large total variation. Con¬ 
sider the general system of conservation laws 

(8) dtU + d,[F{U)]=0, 

where the unknown U{t,x) attains values in , the variables {t,x) G [0,-|-oo[xK 
and the flux function F : —>■ is smooth and strictly hyperbolic (5). Consider 

furthermore the Cauchy problem obtained by coupling ( 8 ) with the initial condition 

(9) U{0,-) = Uo. 

Under some further technical assumption on the structure of the flux, Glimm [16] 
established existence of a global in time solution of the Cauchy problem provided 
that TotVart/o, the total variation of the initial datum, is sufficiently small. Under 
the same assumptions, Bressan and several collaborators established uniqueness 
results, see [5] for a detailed exposition. 

The requirement that the total variation TotVar Uq is small is highly restrictive, 
but necessary to obtain well-posedness results unless additional assumptions are 
imposed on the flux function F. Indeed, explicit examples have been constructed 
of systems and data Uq where TotVar Uq is Hnite, but large, and the admissible 
solution blows up in hnite time. In particular, in [3] Baiti and Jenssen constructed 
an initial datum for system ( 1 ) such that the L°°-norm of the admissible solution 
blows up in hnite time. The solution is admissible in the sense that it is piecewise 
constant and every shock is Lax admissible. For further examples of hnite time 
blow up, see the references in [3] and [14]. 
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2.1.2. Schaeffer’s Regularity Theorem does not extend to systems. In [24] Schaeffer 
established a regularity result which can be loosely speaking formulated as follows. 
Consider a scalar conservation law with strictly convex flux, namely equation (8) 
in the case when U{t, x) attains real values and F : K —)■ K is uniformly convex, i.e. 
F" > c > 0 for some constant c > 0. The work by Kruzkov [20] establishes existence 
and uniqueness of the so-called entropy admissible solution of the Cauchy problem 
posed by coupling (8) and (9). It is known that, even if Uq is smooth, the entropy 
admissible solution can develop shocks, namely discontinuities that propagate in the 
(t, a;)-plane. Schaeffer’s Theorem states that, for a generic smooth initial datum, 
the number of shocks of the entropy admissible solution is locally finite. The word 
“generic” is here to be interpreted in a suitable technical sense, which is related to 
the Baire Category Theorem, see [24] for the precise statement. 

In [II] we discuss whether or not Schaeffer’s Theorem extends to systems of 
conservation laws where every characteristic field is genuinely nonlinear, namely (6) 
holds. Note that the assumption that every characteristic field is genuinely nonlinear 
can be loosely speaking regarded as the analogous for systems of the condition 
(which applies to scalar equations) that the flux is strictly convex. Indeed, regularity 
results for scalar equations with strictly convex fluxes have been extended to systems 
where every characteristic field is genuinely nonlinear: as we mentioned before, this 
is the case of Oleinik’s [23] decay estimate, see for instance [7, 8, 9, 12, 17, 22] for 
possible extensions to systems. Also, the SBV regularity result by Ambrosio and 
De Lellis [1], which applies to scalar conservation laws with strictly convex fluxes, 
has been extended to systems where every characteristic field is genuinely nonlinear, 
see [2, 4, 13]. 

Despite the above considerations, in [11] we exhibit an explicit example which 
rules out the possibility of extending Schaeffer’s Theorem to systems of conservation 
laws where every characteristic field is genuinely nonlinear. More precisely, we 
construct a “big” set of initial data such that the corresponding solutions of the 
Cauchy problems for the Baiti-Jenssen system (1) develop infinitely many shocks 
on a given compact set of the {t, a;)-plane. The term “big” is to be again interpreted 
in a suitable technical sense, which is related to the Baire Category Theorem, see [11] 
for the technical details. 

2.2. The Lax solution of the Riemann problem. We consider a system of 
conservation laws (8) and we assume that F : —>■ is strictly hyperbolic (5) and 
that every characteristic field is genuinely nonlinear, namely (6) holds for i = 1, 2,3. 
Lemma 2.1 below states that the Baiti-Jenssen system satisfies these conditions. 
The Riemann problem is posed by coupling (8) with an initial datum in the form 

( 10 ) 

where U'^ and U~ are given states in K^. In [21], Lax constructed a solution of the 
Riemann problem (8)-(10) under the assumptions that the states [/+ and U~ are 
sufficiently close: we now briefly recall the key steps of the analysis in [21]. 

We fix i = 1, 2, 3 and U G and we define the i-wave fan curve through U by 
setting 

A[s,t/]:=| 


( 11 ) 


i?i[5) fl] S > 0 
F^s,!/] s<0 
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In the previous expression, Ri is the i-rarefaction curve through U and Si is the 
i-Hugoniot locus through U. The i- rarefaction curve Ri is the integral curve of the 
vector field namely the solution of the Cauchy problem 


( 12 ) 


dRi 

ds 


= n{Ri) 


t7] = t/. 


The i-th Hugoniot locus Si is the set of states that can be joined to 17 by a shock 
with speed close to \i{U). The i-Hugoniot locus Si is determined by imposing 
the Rankine-Hugoniot conditions. We term the value |si| strength of the f-wave 
connecting the states U (on the left) and Di[s, U] (on the right). Note that, owing 
to ( 11 ), when Si > 0 the i-wave is a i-th rarefaction wave, when Si < 0 the f-wave is 
an f-shock satisfying the so-called Lax admissibility criterion. The solution of the 
Riemann problem (8)-(10) is computed by imposing 


t/+ =£>3 


S3,D2\s2,Di[si,U ]] 


and by using the Local Invertibility Theorem to solve for (si, 52 , 83 ). From the 
value of (si, S2, S3) one can reconstruct a solution of the Riemann problem (8)-(10), 
see [21] for the precise construction. This solution is obtained by patching together 
rarefaction waves and shocks that satisfy the Lax admissibility criterion. In the 
following, we refer to this solution as the Lax solution of the Riemann problem ( 8 )- 
( 10 ). 


2.3. The wave fan curves of the Baiti-Jenssen system. We collect in this 
paragraph some features of the Baiti-Jenssen system. For the proof, we refer to [3, 
11 ]. 

The first result states that in the unit ball the Baiti-Jenssen system is strictly 
hyperbolic whenever 0 < j] < 1/4. Also, when r; > 0 all the characteristic fields 
are genuinely nonlinear. Note that when rj — 0 this last condition is lost because 
two characteristic fields became linearly degenerate. See [3] or [11] for the explicit 
computations. 

Lemma 2 . 1 . Assume that 0 < rj < 1/4 and that U varies in the unit ball, \U\ < 1. 
Then the Baiti-Jenssen system with flux (2) is strictly hyperbolic, namely (5) holds 
true. If we also have 77 > 0 then every characteristic field is genuinely nonlinear, 
namely (6) is satisfied for i = 1,2, 3. 

We now discuss the structure of the wave fan curves. We start by giving the 
explicit expression of the 1- and the 3-wave fan curve. In the statement of the 
following result, we denote by {u,v,w) the components of the state U G 

Lemma 2.2. Consider the flux function (2), assume that 0 < 77 < 1/4 and fix 
U such that \U\ < 1. Then the following properties hold true. 

i) The Twave fan curve Dfia, U\ is a straight line in the plane v = v, more 
precisely 

(13) D^[a,U] = U + an{U), 

where ri{U) = 0 • 


V 
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Note thatri(U) is the first eigenvector of the Jacobian matrix DF{U). Also, 
the states U (on the left) and Di[a, U\ (on the right) are connected by a wave 
which is 

- a 1-rarefaction wave when cr > 0 , 

— a Lax admissible 1-shock when a < 0. 

a) The 3-wave fan curve D^ItjU] is a straight line in the plane v = v, more 
precisely 


(14) 


D^[t, U] = U + tUU), 


where r^iU) = 



The vector r^(U) is the third eigenvector of the Jacobian matrix DF{U). 
Also, the states U (on the left) and D^It, U] (on the right) are connected by 
a wave which is 

- a 3-rarefaction wave when r < 0, 

- a Lax admissible 3-shock when t > 0. 


Note that, for the 3-wave fan curve, the positive values of r correspond to shocks, 
the negative values to rarefaction waves. This is the contrary with respect to (11) 
and it is a consequence of the fact that we use the same notation as in [3, 11] and we 
choose the orientation of r^ in such a way that when 77 > 0 condition ( 6 ) is replaced 
by the opposite inequality 

VA 3 • ra < 0. 

We now turn to the structure of the 2-wave fan curve. In the following statement, 
we use the notation 



Also, we consider entropy admissible solutions of scalar conservation laws, in the 
Kruzkov [20] sense. 


Lemma 2.3. Assume that U is a Lax solution of the Riemann problem ( 8 )-(10). 
Then the second component v is an entropy admissible solution of the Cauchy prob¬ 
lem 


(15) 


dtv -\- = 0 

^'(O^a;) = I 


X < 0 
X > 0. 


Also, we can choose the eigenvector r 2 and the parametrization of the 2-wave fan 
curve D 2 [s, U] in such a way that the second component of D 2 [s, U] is exactly u -I- s. 


3. Interaction of two 2-shocks 


We first rigorously state Lemma 1.1 

Lemma 3.1. There is a sufficiently small constant e > 0 such that the following 
holds. Fix a constant a such that 0 < a < 1/2 and set C/jj := (a, 0, —a). Assume 
that 

\Ue — U^\ < ea, 0 < p < ea, 

Si,S2 < 0, Si,S2 e [-£a,0]. 
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Assume furthermore that 


(16) 


Ur = D2 S2,-C>2[si, C4] 


Then there are a < 0 and r > 0 such that 


(17) 


Ur=D3 T,D2 [si + S2,Di[a,Ug]] 


Note that by combining (17) with the inequalities cr < 0, r > 0 and si + S 2 < 0 
we get that the three outgoing waves are all shocks. The proof of Lemma 3.1 is 
organized as follows: 

§ 3.1: by relying on a perturbation argument, we show that the proof of Lemma 3.1 
boils down to the proof of the Taylor expansion (21). 

§ 3.2: we complete the proof by establishing (21). 

3.1. Proof of Lemma 3.1: first step. We start with some preliminary consid¬ 
erations. Assume that the states Ug and Ur satisfy (16). Next, solve the Riemann 
problem between Ug (on the left) and Ur (on the right): owing to [21], this amounts 
to determine by relying on the Local Invertibility Theorem the real numbers tr, s 
and T such that 


Ur=D3 T,D2[s,Di[(T,Ug]] 


(18) 


Establishing the proof of Lemma 3.1 amounts to prove that s = si -I- S 2 <0 and 
that (T < 0 , r > 0 . 

To prove that s = si -I-S2 we recall Lemma 2.3 and the fact that the v component 
is constant along the 1-st and the 3-rd wave fan curves Dg and D3. We conclude 
that s = — Uf = Si -f S 2 < 0 . Note that Vr and vg are the second component of 

Ur and Ug. 

We are left to prove that cr < 0 and t > 0. We first introduce some notation: we 
regard tr and r as functions of rj, si and S 2 and Ug and we write cr^(si, S 2 , Ug) and 
Triisi, S 2 ,Ug) to express this dependence. Note that a and r depend on ij because 
the wave fan curve D2 depends on ij. 

Owing to the Implicit Function Theorem, the regularity of cr^(si, S 2 , t^) and 
o’riisi, S2,Ug) is at least the same as the regularity of the functions Dg, D2 and 
D 3 . Also, note that the Lax Theorem [21] (see also [5, p.lOl]) states that the wave 
fan curves Dg, D2 and D3 are C^. The reason why we can achieve C°° regularity 
is because we are actually considering the wave fan curves in regions where they 
are C°°. To see this, we first point out that, owing to (13) and (14), the wave fan 
curves Di, are straight lines and hence they are C°°. Next, we point out that we 
are only interested in negative values of si -I- S 2 . Hence, we can replace the 2-wave 
fan curve D2 defined as in (11) with the 2-Hugoniot locus 82- We recall that the 
2-Hugoniot locus S' 2 [s,C/] contains all the states that can be connected to ?7 by a 
shock, namely all the states such that the couple {U, S' 2 [s, C7]) satisfies the Rankine- 
Hugoniot conditions. The 2-Hugoniot locus S' 2 [s, 17] is C°° and by combining all 
the previous observations we can conclude that <Jn{si,S2,Ug) and S2,Ui) are 
both C°° with respect to the variables (ry, si, S 2 , Ug). 

Next, we discuss the partial derivatives of crjj(si, S2,Ug) and Tjj{si, S2,Ug) with 
respect to (si,S 2 ) at the point ( 77 , 6 , 0 , lA). By arguing as in the proof of estimate 
(7.32) in [5, p.l33] we conclude that 
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(19) 


for every Ut, for every 77 > 0 and every integer A: > 1 we have the following 
equalities: 



d^ar, 

d'^Trf 

d’^Tr, 

ds’i 

(0,0,Ut) 

(0.0.;7f) 

(0,0,Ut) ^'^2 


For every Ui and for every 77 > 0 we also have the following equality con¬ 
cerning the derivatives of second order: 


dsids 2 


( 0 , 0 ,Ue) 


d'^Tn 

dsids 2 


(0,0,Ut) 


= 0 . 


This implies that and admit the following Taylor expansions 


( 20 ) 


C’'r)(Si, S 2 , Ul) 


(^1: "^2; ) 


1 9 V^( 0 , 0 , 17 f) 2 , i a%( 0 , 0 ,f/f) 2 

2 ds\ds2 2 dsids\ 

-|-0(|(S1,S2)|)S1S2(S1 -I-S2) 

15 %( 0 , 0 ,C/^) 2 , ia%(o,o,c/,) 2, 

2 " d4ds2 . + 2 " dsM . 

-|-o(|(si,S2)|)siS2(si +S2) 


In § 3.2 we prove that when 77 = 0 and Ui = the functions cr and r admit the 
Taylor expansions 

( 21 a) ^ (_\)«i« 2 (si +S 2 ) +o(|(si,S 2 )n. 

Next, we use the Lipschitz continuous dependence of the derivatives of third order 
with respect to 77 and and we conclude that 

ds\ds2 16 

g^r^( 0 , 0 , Ut) 
ds\ds2 16 

provided that 0 < rj < ea and |[/f — { 7 jj| < ea. In the above expression, C denotes a 
universal constant. By plugging the above expressions into (20) and recalling that 
Si, S 2 < 0 we can eventually conclude that, if e is sufficiently small, then 

ar,{si,S2,Ul) < ^SiS2{si + S2) < 0 , 

64 

'b)(si, S2, Ui) > ——SiS2(si + S2) > 0 . 

64 

The proof of the lemma is complete. 


+ 


a%(o,o,c/,) 


dsids 2 

d^Tn{ 0 , 0 , Ui) 


dsidso 


a 

' 16 
a 
16 


< Cea 


< Cea 


3.2. Proof of formula (21). The proof of the Taylor expansion (21) is divided 
into two parts: 

§ 3.2.1: as a preliminary result we determine the structure of the Hugoniot locus 
^ 2 [S,C/] 

§ 3.2.2: we conclude the proof. 

Note that in this paragraph we always assume 77 = 0 because formula (21) deals 
with this case. 
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3.2.1. The 2-Hugoniot locus. Before giving the technical results, we introduce some 
notation. First, we recall that we term Fq the flux function in (2) in the case 
when r] = 0. In the following, we will mostly focus on the behavior of the first and 
the third component of U. Hence, it is convenient to term U and Fq the vectors 
obtained by erasing the second components of U and Fq, respectively. We have the 
relation 

where we have also introduced the 2x2 matrix J{v). 

Finally, we recall that we term 5 ' 2 [s, U] the 2-Hugoniot locus passing through U, 
namely the set of states that can be connected to C/ by a (possibly not admissible) 
shock of the second family. Also, as usual we denote by u, v and w the first, second 
and third component of t/, respectively. We use the notation U = {u,w). 


Lemma 3.2. Fix rj = 0 and assume that \2v + s| <4, then the 2-Hugoniot locus 
through U has the following expression: the second component of 52 [s, 17] is v s 
while the first and third components are 

(23) S 2 [s,U]=U + Eiv,s)U 


where the 2x2 matrix E(v,s) is 

4s / s + 4-2u 

~ (2u + s)2 - 16 V(s + 4)(s - 2) + 4u 


3s — 4 + 2uy 


Proof. By Lemma 2.3 the second component of 52 [s, C/j is u + s. To construct 
52 [s, U] we use the Rankine-Hugoniot conditions, which are a system of 3 equations. 
Owing to Lemma 2.3, the second equation reads 

7S = (u + s)^ — 

and this implies that the speed 7 of the 2 -shock is 

(24) "f = 2v s. 

We define the vector 2l(s, U) by setting 

2i{s,U) ■■= K[s,U]-U 

and we point out that to establish Lemma 3.2 we are left to show that 

(25) 2i{s,U) = E{v,s)U. 

The first and the third equations in the Rankine-Hugoniot conditions can be written 
as 


(26) 


'y2l{s,U) = J(v + s) U + 2l{s,U) 


J{v)d, 


where J is the same as in (22). Next, we introduce the 2x2 matrix 
A{v, 7 ) = 7 J — J(y) 

^7 0^ 

.0 Ij 


-4 


u — 1 —1 

v{v — 2 ) 1 — v) ' 


and we rewrite (26) as 

A (u -h s, 7 ) 2t(s, U) = [j(u -h s) - J(u)l U, 
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which implies (25) provided that 

E{v, s) = A^^{v + s, 7 ) [J{v + s) — ^(■p)] 

By recalling that 7 = 2v + s we can compute the explicit expression of the above 
matrices: 

A{v + s,2v + s)= (^^^2 _y_ + s) 61; + 5s - 4 ) ’ 

j(^ + ,)_ j(^)=4s(^^^2l-2 -l)- 

The determinant of the matrix A(y + s, 21; + s) is 

det := ( 2 i; + s)^ — 16 

and hence the matrix is invertible when |2i; + s| <4. We can now complete the 
lemma by computing the explicit expression of E, namely 

r(- 6 u + 5 s -4 -4 \ / 1 0 \ 

^y4:{v + s-2){v + s) 4:-3s-2v)^ \s + 2v-2 -ij 

_4s/s + 4-2i; □ 

det \(s + 4)(s — 2) +4v 3s — 4 + 2vJ ' 


3.2.2. Conclusion of the proof of formula (21). We are now ready to establish (21). 
We first recall some notation: we consider the system of conservation laws with 
flux Fq, see (2). We consider the collision between two 2-shocks and we assume 
that t/jj = (a, 0,—a), Um and Ur are the left, middle and right states before the 
interaction. This means that for some Si < 0, S 2 < 0 we have 


(27) 


Ur — D2 [s 2, Um] — D 2 [s 2 , 7l2[si, C^j]] 

^ S 2 [S 2 ,S 2 [SUU{\]. 


In the above expression, S 2 represents the 2-Hugoniot locus. To establish the last 
equality we used the fact that si and S 2 are both negative. We plug (23) into (27) 
and we use the equality = 0 : we arrive at 


(28) 


Ur — Ufj + E[0, si)U^ +E{si,S 2) U\j + E{0, si)U^ 
= U(^ + E{0^ si) + E(si, S 2 ) 4-E[si^ S2)E{0, si) U^ 


Next, we focus on the states after the interaction. By arguing as at the begin¬ 
ning of § 3.1, we conclude that it suffices to determine a = cro(si, S 2 , t/j) and 
T = To(si, S 2 , t/j) such that 


Ur=D3 


t,D2 [si -\- S2,Di[a, Ujj]] 


By the explicit expression of Di and D 3 and by applying Lemma 3.2 we infer that 
the above equality implies 


Ur = 


Ujj-|-crri(O) -f £1(0, Si -|- S2) t/(| -b crri(O) + Tr 3 (si-|- S2) 


(29) 


— 17(1 -b £1(0, Si -b S2)C7(j -b £ -b £1(0, Si -b S2) (Tri(O) -b Tr^{si -b S2) 


— 17(j -b £1(0, Si -b S2)C7(j -b H(si -b S2) 
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In the previous expression we denote by ri and ra the vectors obtained from ri 
and r 2 by erasing the second component. Also, we introduced the matrix H: its 
first column is [/ + £1(0, si + S 2 )]ri( 0 ), the second column is r 3 (si + 52 )- In the 
following, we will prove that H(si + S 2 ) is invertible provided that si and S 2 are 
both sufficiently close to 0. By comparing (28) and (29) we then obtain 
(30) 


— H ^(si + S2) £1(0,Si) + £/(si,S 2 ) + £1(31,S2)£1(0,Si) — £1(0,Si + S 2 ) 


G^(si, S2) 


Assume that we have established the following asymptotic expansion for G: 

(31) ~ ^ 2 3 ^ *1*2(51 + S2) + o(|(si, S2)P). 

Then by plugging both (31) and = (a, —a) into (30) we obtain the asymptotic 
expansion (21). Hence, to conclude the proof of (21) we are left to establish (31). 

First, we point out that, owing to the expression of E in the statement of 
Lemma 3.2, 

(0>*)-g2_i6 (s + 4)(s-2) 3s-4j- 

This implies that when Si = S 2 = 0, the matrix £1(0, Si + S 2 ) vanishes and hence 

H-(0)=(n(0)|,i,(0))"™i'-‘>(; ^ 2 ) ' = (; % )■ 

We compute now the asymptotic expansion of 


-E1(0, si) + £1(0 + si, S 2 ) — £1(0, Si + S 2 ) + £1(0 + si, S2)£l(0, si). 


By directly computing the sum of the above matrices, we obtain that we can factor 
the term 

_ 4 siS 2 (si + S2)_ 

(sf - 16 )((si + S2)2 - 16 )(( 2 si + S2)2 - 16 ) ’ 
which multiplies the matrix with coefficients 

Coeffi,i: (si + 4 )(si + S2 + 4 )( 6 si + 5 s 2 — 12 ) 

Coeffi^2: 4(5 s^ + 13siS2 + 9s^-48) 

Coeff2T: 2(si + 4)(si + S2 + 4)(4 - 6si 

+ 2si — 7 s2 + 4siS2 + 2 S 2 ) 

Coeff2.2: 192 - 128si - 36s? + 26s? 

— 160 s 2 — 52siS2 + 65s?S2 — 20s? 

+ 55sis? + 16 s?. 

By combining the above computations we obtain the following asymptotic expan¬ 
sion: 

G'(si, S 2 ) =^ Q -1/2 ) ( 2-4^ 3-4^ ^ • SiS2(si + S 2 ) + o(||(si,S2)P) 

( 32 ) =4(2 3 )*l*2(*l+*2)+0(||(si,S2)f). 
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This establishes (31) and hence concludes the proof of (21). 


□ 


4. Interaction of a 1-shock and a 2-shock 


We first rigorously state Lemma 1.2. 


Lemma 4.1. There is a sufficiently small constant e > 0 such that if 0 < rj < e, 
then the following holds. Assume that the states Ui, Ur G satisfy 


(33) 


Ur = Di 


a, D 2 [s,Ue] 


for real numbers s, a such that 

(j,s < 0, |s|, \a\ < 

Furthermore, assume that \Ui\ < 1/2. Then there are real numbers a' and t' such 
that 


(34) 
and 

(35) 


Ur = D. 


t' ,D2 [s,Di[a' ,Ui]] 


2a < a' < -cr, ^ lOcrs. 


Note that (35) implies a' < 0 and r' > 0. If we combine these inequalities 
with (34) and s < 0 we see that the three outgoing waves are all shocks. 
Establishing the proof of Lemma 4.1 amounts to establish (35). Indeed, 

(1) by using Lax’s construction (see § 2.2) we determine a' ,s’ ,t' such that 


Ur = D^ 


r',L >2 [s ,Di[a', Ui]] 


(2) By combining (13), (14) and Lemma 2.3 we obtain that s' = s. 

To establish (35) we proceed as follows: 

§ 4.1: we establish (35) in the case when 77 = 0. 

§ 4.2: we conclude the proof by relying on a perturbation argument. More pre¬ 
cisely, by using the fact that the flux E/ in ( 2 ) smoothly depends on rj we 
show that (35) holds provided t] is sufficiently small. 

Note that, as we have mentioned in the introduction, the precise expression of the 
function pi and ps plays no role in the proof of Lemma 4.1, what is actually relevant 
is that estimates (35) hold at 77 = 0 with strict inequalities and that 77 is sufficiently 
small. 


4.1. Proof of Lemma 4.1: the case 77 = 0. We establish (35) in the case rj = 0. 
This part of the proof is actually the same as in [3, p. 844-845], but for completeness 
we go over the main steps. 

We term Cg, r/ the real numbers satisfying (34) when 77 = 0. Let Vr and V£ 
denote the second components of Ur and U^, respectively. We term 7 the speed 
of the incoming 2 -shock (which is the same as the speed of the outgoing 2 -shock), 
we recall Lemma 2.3 and the fact that the second component varies only across 
2-shocks. We conclude that 

9 2 

Vt- Vp 

7 = - = Vr + Vl = 2V£ + S. 

Vr - Vp 

Since by assumption jt/^j < 1/2 and jsj < 1/4, then 
(36) |7l < 3. 
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By imposing the Rankine-Hugoniot conditions on the incoming and outgoing 2- 
shocks and by arguing as in [3, pp. 844-845], with the choice c = 4, we arrive at 
the following system: 

(7 + 4)a'o + (j- 4)t^ = (7 + 4)cr 
vei'y + 4)(Tq -I- s - 2)(7 - 4)rg = {vi -|- s)(7 -I- 4)(t 


If we set 

(37) 
and 

(38) 


A := 


Xn = 


7-1-4 7 — 4 

V£(j + 4) (v£ + s - 2)(j - 4) 


r = 


(ve + s){'y + 4) 


then the above linear system can be recast as AXq = Y. The explicit expression of 
the matrix A~^ is 


(39) 


1 


(40) 


(42 _72)(_ s -^2) 

wi 
2 


{v£ + s - 2){-f - 4) -( 7 - 4 ) 

-z;^(7-|-4) (7 + 4 ) 


We solve for dg and Tq and we obtain 


—s -t- 2 


Tn = 


(4-7)(-s + 2) 


sa 


By using (36) and the inequality jsj < 1/4, we obtain 


2 
- < 


1 

< 1 , XT < 


7- 


< 4 


3"-s-b2""’ 21 " ( 4 - 7 )(-s-k 2 ) 

and this implies that the estimate (35) holds true in the case when 77 = 0. 


4.2. Proof of Lemma 4.1: the case rj > 0. We are now ready to complete the 
proof of Lemma 4.1. We proceed as follows: 

§ 4.2.1: we make some preliminary considerations which reduce the proof of Lemma 4.1 
to the proof of the fact that a certain map is a strict contraction. 

§ 4.2.2: we conclude the proof by showing that the map is indeed a strict contrac¬ 
tion. 


4.2.1. Preliminary considerations. We first introduce some notation. We term Um 
the intermediate state before the interaction, namely 

(42) := D2[s,Ue]. 

Also, we term and the intermediate states after the interaction, namely 

U^:=D,[a',Ue], 

(43) 17" := [s, U'J = ila [ - r', 1 / 3 ] = [ - r', 111 [a, Um]] 

^D3[-T',D,[a,D2[s,Ue]]] 

Next, we use [3, eq. (5.3)-(5.4)] and we recast the Rankine-Hugoniot conditions for 
2 -shocks as a nonlinear system in the form 

(44) AX + r]T{X,Ue,s,a) = Y, 

where A and Y are as in (37) and (38), respectively. Also, the vector X is defined 
by setting 


a: := 


r 
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and the nonlinear term T{X, Uf, s, a) is equal to 


{ Pl{U^^)-Pl{Um)-pi{Ui^)+Pl{Ul) \ 

{p3(.Ui^)-p3iUm)-p3iU^^)+P3m, )■ 


In the above expression, the functions pi and ps are the same as in (3). Note, 
however, that the precise expression of pi and pa plays no role in the proof, the 
only relevant point is that pi and pa are both regular (say twice differentiable with 
Lipschitz continuous second derivatives). Note furthermore that we can regard X 
as a function of X, Ug, s and tr because, owing to (42) and (43), Um, U!^ and U'^ 
are functions oi X, Ui, s and cr. Next, we rewrite equation (44) as 

(46) X = Xo - vA -^ XiX , Ui , s , a ), 

where the vector Xq = A~^Y is given by (38) and (40). 

We now fix s, a, p and \Ui\ satisfying the assumptions of Lemma 4.1 and we 
define the closed ball 


(47) A:= {X = (cr',/) e . |X-Xo| < kpas}. 

In the above expression. A: > 0 is a universal constant that will be determined 
in the following and ctq and Tq are defined by (40). We also define the function 
T : by setting 

(48) TiX):=Xo-pA-^XiX,Ui,s,a). 

Assume that T is a strict contraction from A to A. Then the proof of Lemma 4.1 is 
complete: indeed, owing to (46) the fixed point X satishes the Rankine-Hugoniot 
conditions (44). Also, owing to (41) and to (47) we infer that the inequalities (35) 
are satisfied provided that the parameter p is sufficiently small. 


4.2.2. Conclusion of the proof of Lemma 4-1. In this paragraph we prove that the 
map T defined by (48) is a strict contraction on the closed set A defined by (47). 

First, we make some remarks about notation. To simplify the exposition, in 
the following we denote by C a universal constant: its precise value can vary from 
occurrence to occurrence. Also, in the following we will determine the constant k 
in (47) and then choose the constant p in such a way that kp < 1. This choice 
implies in particular that, when X belongs to the set A defined as in (47) and the 
hypotheses of Lemma 4.1 are satisfied, then the map X attains values on a bounded 
set and so X and all its derivatives are bounded by some constant C. Finally, note 
that, X & A and the hypotheses of Lemma 4.1 are satisfied, then |A“^| < C. 

We now proceed according to the following steps. 

Step 1: we point out that to show that the map T is a contraction it suffices to 
show that 

(49) \X{X,Ut,s,a)\<Cas 

provided that X & A and the hypotheses of Lemma 4.1 hold. Indeed, assume 
that (49) holds, then 

(48) (49) 

\T{X)-Xo\ < p\A-^F{X,Ui,s,a)\ < Cpas 
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and hence T attains values in the set ^ defined as in (47) provided that k is large 
enough. Also, 

(48) 

|T(Ai)-T(A 2 )| < r,\A-^\\J^{Xi,Ue,s,a) - J^iX 2 ,Ue,s,a)\ 

< riC\Xi — A2I < 2 1^1 ~ ^2! 

provided that the constant 77 is sufficiently small. This implies that T is a contraction 
and concludes the proof of Lemma 4.1. 

Step 2: we establish (49). First, we point out that, if A € A and the hypotheses 
of Lemma 4.1 are satisfied, then 

bi(C/") -Pi(t7™) -^1(1/;,) +Pim\ < c(|t/" - u^\ + \u:^ - Ue\) 

< c[\u';,-uA + \Ur-u^\ + \uU-u,\) 

< c{\t'\ + |a| + |a'|) < c{\a\ + |A'| + |A - A'|) 

(40). (47) / s 

< C\^\a\ + rjkasj < C|cr|. 

By using an analogous argument, we control the second component of F and we 
arrive at 

(50) \F{X,Uus,a)\<C\a\. 

Next, we point out that when s = 0 we have U'^ = and Ui = Um and by using 
again the Lipschitz continuity of the functions pi and p^ we conclude that 

(51) |A(A,C/^,s,a)| <C|s|. 

Finally, we use the regularity of the function F and, by arguing as in the proof of [5, 
Lemma 2.5, p. 28], we combine (50) and (51) to obtain (49). This concludes the 
proof of Lemma 4.1. □ 
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